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Abstract 

In this paper diffusion processes with changing modes are studied involving the 
variable order partial differential equations. We prove the existence and uniqueness 
theorem of a solution of the Cauchy problem for fractional variable order (with 
respect to the time derivative) pseudo-differential equations. Depending on the 
parameters of variable order derivatives short or long range memories may appear 
when diffusion modes change. These memory effects are classified and studied 
in detail. Processes that have distinctive regimes of different types of diffusion 
depending on time are ubiquitous in the nature. Examples include diffusion in a 
heterogeneous media and protein movement in cell biology. 
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1 Introduction 

Diffusion processes can be classified according to the asymptotic behaviour of their mean 
square displacement (MSD) as a function of time. If the dependence of the MSD on 
time is linear then the process is classified as normal, otherwise as anomalous. For many 
processes, the MSD satisfies 

MSD(t) ~ Kpt 13 , t^oo, (1) 

where Kp is a constant. If /3 = 1 the diffusion is normal, if (3 > 1 the process is super- 
diffusive, while if (3 < 1 the process is sub-diffusive [2TJ [36]. The ultra-slow diffusion 
processes studied in [5j [20j 122] lead to logarithmic behaviour of MSD for large t. The 
MSD of more complex processes with retardation (see [51 [12]) behaves like t^ 2 for t 
small, and t^ 1 for t large, where (3% < /3 2 . Subdiffusive motion with 0.1 < (3 < 0.9 was 
recorded in [TO], and with 0.22 < (3 < 0.48 in [29], depending on macromolecules and 
cells. In [TBI [T7] protein movement is studied in the cell membrane with a few types of 
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compartments and made a conclusion that (3 depends on time scales. Our models are 
subdiffusive, but of variable order with order function (3 = (3{t). 

It is well known that simple homogeneous subdiffusive processes can be modeled 
using a fractional order partial differential equation where only the time derivative has 
a constant fractional order [2T]. Variable fractional order derivatives and operators 
were studied by N. Jacob, at al. [15], S. G. Samko, at al. [23 EH], W. Hoh [H]. 
Recently A. V. Chechkin, at al. [6] used a version of variable order derivatives to 
describe kinetic diffusion in heterogeneous media. In the recent paper [19], Lorenzo 
and Hartley introduced several types of fractional variable order derivatives and applied 
them to engineering problems. We will modify these operators, restrict them to order 
functions (3{t) that are piecewise constant and then apply the resulting variable order 
partial differential equations (VOPDE) to diffusion processes with changing diffusion 
modes. An important aspect of the modeling is that the VOPDEs provides a description 
of memory effects arising from a change of diffusion modes that are distinct from the 
"long range memory" connected with the non-Markovian character of diffusion. Thus, 
in the VOPDE based description of anomalous diffusion models, both non-Markovian 
long range memory and new type of memory may be present simultaneously. 

The paper is organized as follows. In Section 2 we introduce background material. In 
Section 3 we study the memory effects arising in connection with a change of diffusion 
modes. In Sections 4 and 5 we study the mathematical model of diffusion processes 
with changing modes in terms of an initial value problem for VOPDE. Namely, we prove 
the theorem on the existence and uniqueness of a solution of the initial value problems 
for variable order differential equations and study some properties of a solution. The 
theorems are proved under the assumption that diffusion mode change times are known. 



2 LH-parallelogram and variable order derivatives 

Recently Lorenzo and Hartley [19j introduced three types of derivatives of variable frac- 
tional order (3(t), t > 0, < (3(t) < 1 all of which are special case of a more general 
fractional order derivative 

&$m = ^f Q &${t,T)f{T)dT, (2) 

where \i and v are real parameters, t > 0, and 

/Cf ( * ) (t,r) = — -sp - r , 0<r<t. (3) 

For convenience in studying of initial value problems, we prefer to use the closely related 
Caputo type operator 

f *®&T)^dT. (4) 
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To describe the properties of the kernel ([3]) and the fractional derivative operators 
(|2j) and (jlj) we introduce the Lorenzo- Hart ley (LH) causality parallelogram [19] II = 
{(/x, v) e R 2 : < fi < 1, -1 < i/ < +1, < fi + v < 1} . The kernel ©, and thus, both 
the operators (J2J) and (jlj) are weakly singular for (//, v) G II. Further, denote 

/C ( t ' r ' s ) = ra-^Kt-r)^) ' * > °' < r < *' s ~ °' (5) 

where < (3(s) < 1 Q is a given function, which is called an order function. 

Our main goal is to model problems where for different time intervals there are 
different modesj of diffusion. To this end, let Tj be a partition of the interval (0, oo) 
into N + 1 sub-intervals (Tk,Tk+i), where = To < T\ < ... < Tjv < Tjv+i = oo. Then 
let P(t) be the piecewise constant function 

N 

P(t)= y %2/3 k I k (t), te(0,oo), (6) 

k=0 

where X k is the indicator of the interval (Tjt,Tt + i) and < l3 k < 1, k = 0,...,N, are 
constants. Under these conditions, the function (j5J) becomes 

N 1 

K(t, r, S ) = Y,Us) T{1 _ m _ Trk , t > 0, < r < t, s > 0, (7) 

and the kernel of the fractional order operator (jlj) becomes 

/C2J } (t, r) = #(t, r, /it + vt), t > 0, < r < t (8) 

with K(t,r,s) defined in (J7j). 

We think of the input to our model as the triplet is), < k < N, while the 

output of our model is determined by the kernel (jSj). Correspondingly, we say that 
the triplet [f3 k ,ii, u) determines the diffusion mode in the time interval (T k ,T k+ i). The 
output is determined by which values of (3{t) are used to compute the variable order 
derivative, that is, by which interval (T k , T k +i) the point fit + vt belongs to. We always 
assume that (/i, v) £ II and then note that r e (0, t) yields fit + vt 6 (fit, (fi + v)t) 
and that (fit,(fi + v)t) C (0,i). This means that the operators T>^§ and Pf^*], use 
information taken in the time sub-interval (fit, (fi + v)t) if v is positive and from the 
sub-interval ((fi + v)t,fit) if v is negative. In both cases, the length of this interval is 
\v\ t. The condition (fi, v) e II predetermines the causality, since < fit + vt < t for all 
t > and < r < t. 

2.1 Generalized function spaces ^ Gp (iR 71 ), ^'_ Gq {R n ) 

Let p > 1, g > 1, p -1 + = 1 be two conjugate numbers. The generalized functions 
space ^_ Gq (]R n ), which we are going to introduce is distinct from the classical spaces 

Mf p(t) = 1, then we agree V^%f(t) = 
2 see definition in Section [2731 
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of generalized functions. In the particular case of p = 2 this space was first used by 
Yu.A.Dubinskii [8j in the course of initial-value problems for pseudo-differential equa- 
tions with analytic symbols. Later, the general case for Vp was studied in [311 [32]. Here 
we briefly recall some basic facts related to these spaces, referring the interested reader 
to [3H 03] for details. 

Let G C M n be an open domain and a system Q = {gk}kLo of open sets be a locally 
finite covering of G, i.e., G = {J^L gk, 9k CC G. This means that any compact set 
K C G has a nonempty intersection with a finite number of sets gu- Denote by {4>k}^L 
a smooth partition of unity for G. We set Gn = Uf? = i9k an d «jv(£) = Y^k=\ ■ ^ 
is clear that Gn C Gn+i, N = 1, 2, ... , and Gn — » G for iV — > oo. Further, by F/ 
(or /(£)) for a given function /(x) we denote its Fourier transform, formally setting 
Ff(0 = Jut™ f(x)e lx ^dx, and by F^ 1 f the inverse Fourier transform, i.e. F~ l f(x) = 
(2n)~ n J Kn f(£)e~ tx ^d£. The support of a given / we denote by supp/. 

Definition 2.1 Let N G W. Denote by ^n,p the set of functions f G L p (M n ) satisfying 
the conditions (1)—(S): 

(1) supp Ff c G N ; 

(2) supp Ff fl supp <pj = for j > N; 

(3) PnU) = WF^KNFfWp < oo. 
Lemma 2.1 For N = 1, 2, ... , the relations 

^n, p ^ ^n+i, p , ^n, p ^ L p (R n ) 
hold, where denote the operation of continuous embedding. 

It follows from Lemma 12.11 that \I/ n,p form an increasing sequence of Banach spaces. 
Its limit with the inductive topology we denote by *&g,p- 

Definition 2.2 ty GtP (M n ) = ind lim^oo ty NtP . 

The inductive limit topology of ^G tP (M n ) is equivalent to the following convergence. 

Definition 2.3 A sequence of functions f m G \£ , G iP (iR n ) is said to converge to an ele- 
ment f G M/ ( ;,,( //,'") iff: 

1. there exists a compact set K C G such that supp f m C K for all m G IN ; 

2- \\fm ~ fo\\ P = (J Rn \ fm ~ fo\ p dx)p > for m > oo . 

Remark 2.1 According to the Paley-Wiener-Schwartz theorem, elements of ^ c,p{JR n ) 
are entire functions of exponential type which, restricted to M n , are in the space L p (M n ). 

The space topologically dual to \1/G,p(-K n ), which is the projective limit of the se- 
quence of spaces conjugate to ^n, p , is denoted by \I/_ G q (M n ). 
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Definition 2.4 ^'_ Gq {lR n ) = pr lim^*, ^ 

In other words, ^>'_ Gq {IR n ) is the space of all linear bounded functionals defined on 
the space ^a,p(Fi n ) endowed with the weak topology. Namely, a sequence of generalized 
functions G ty_ G (M n ) converges to an element g G ty_ Gq (M n ) in the weak sense 
if for all / G ^a,p(^R n ) the sequence of numbers < gw, f > converges to < g , f > as 
N — > oo. We recall that the notation < g, f > means the value of g G ^ _ G q (M n ) on an 
element / G ^ GjP (iR n ). 

2.2 Pseudo-differential operators with constant symbols. 

Now we recall some properties of pseudo-differential operators with symbols defined and 
continuous in a domain G C M n . Outside of G or on its boundary the symbol A(£) may 
have singularities of arbitrary type. For a function <p G ^G,p(^ n ) the operator A(D) 
corresponding to A(£) is defined by the formula 

A(DMx) = j^y n j G A{Z)F<p{Z)e**dt. (9) 

Generally speaking, A(D) does not make sense even for functions in the space C^{IR n ). 
In fact, let £o be a non-integrable singular point of A(£) and denote by O(£o) some 
neighborhood of £ Q . Let us take a function ip G C™(M n ) with F(p(£) > for £ G O(£ ) 
and F<p(£o) = 1- Then it is easy to verify that A(D)ip(x) = oo. On the other hand, 
for <p G \&G,p(-K n ) the integral in Eq. (jHJ) is convergent due to the compactness of 
supp Fip C G. We define the operator A(—D) acting in the space ^_ Gq (M n ) by the 
duality formula 

< A(-D)f, <p> = <f, A(D)<p >, / G *'_ Gi9 (i?n ip G * Gj ,(^). (10) 



Lemma 2.2 T7ie spaces ^/c,p(^i n ) and ty'_ G JlR n ) are invariant with respect to the ac- 
tion of an arbitrary pseudo- differential operator A(D) whose symbol is continuous in G. 
Moreover, if A(£)kn(£) is a multiplier on L p for every N G IV, then this operator acts 
continuously. 

Remark 2.2 In the case p = 2 an arbitrary pseudo-differential operator whose symbol 
is continuous in G acts continuously without the additional condition for A(£)kjv(0 t° 
be multiplier in L 2 for every N G IN. 

2.3 Subdiffusion processes. 

As is known [HI EI], a (sub-) diffusion process is governed by the fractional order partial 
differential equation 

D%u(t, x) = A(D)u(t, x), t > 0, x G R n , (11) 
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where D% is the Caputo fractional derivative of order (3 G (0,1], and A(D), D = 
(Di, ...,D n ), Dj = —i-£r, j = 1, ...,n, some elliptic pseudo-differential operator. 

Many diffusion processes driven by a Brownian motion can be described by equation 
(fTTT) with a second order elliptic differential operator A(D) and (3 = 1. Levy stochas- 
tic processes (which include jumps) also connected with (TTTj) and an elliptic pseudo- 
differential operator A(D) (see, e.g. [2]). In particular, if particle jumps are given by a 
symmetric Levy stable distribution with infinite mean square displacement then A(D) 
is a hyper-singular integral, defined as the inverse to the Riesz-Feller fractional order 
(0 < a < 2) operator (for details see [26J ) . A wide variety of non-Gaussian stochastic 
(subdiffusive) processes lead to equation ffTTj) with < (3 < 1 (see [2U [22]). For diffu- 
sion governed by distributed order differential equations see [H [31] . The parameter f3 
determines the sub-diffusive mode, which is slower than the classical free diffusion. 

Generalizing this approach we will say that the fi, z/}-diffusion mode in the time 
interval {Tk-i, Tjt) is governed by the equation 

V^ u} u{t, x) = A{D)u{t, x),te (T fc _!, Tk), x G M n . (12) 

The entire process then can be described by the equation 

V *{lu} u (^ x ) = A(D)u(t, x),t>0,xeM n , (13) 

where T)^j£ \ v \ is the variable fractional order operator with the kernel K?^K(t, r) in (jSJ). 

In Section 4 we will prove the existence of solutions to the initial value problem 
defined by the differential equation ffl3]) . 

3 Changing of modes: 'short-range' and 'long-range' 
memories 

We call the triplet {Pk,[i,u} admissible if < (3f. < 1 and (/i, v) G II. Diffusion in 
complex heterogeneous media is accompanied by frequent changes of diffusion modes. 
It is known that a particle undergoing non-Markovian movement possesses a memory 
of past (see [211 [36]). Protein diffusion in cell membrane, as is recorded in [211 [25] is 
anomalous diffusion. Descriptions of this process using random walks also shows the 
presence of non-Markovian type memory [UdSllIS]- It turns out, there is another type 
of memory noticed first by Lorenzo and Hartley in their paper [19] in some particular 
cases of // and v. This kind of memory arises when the diffusion mode changes. 

In this section we study a special case of this phenomenon where there is a single 
change of diffusion mode, that is, a sub-diffusion mode given by an admissible triplet 
{Pi,/jl,v} changes to a sub-diffusion mode corresponding to another admissible triplet 
{/3 2 , fi, u} at some particular time T. 

Definition 3.1 Let v} and {02, be two admissible triplets. Assume the 

diffusion mode is changed at time t = T from {Pi, fi, v}-mode to {/3 2 ,/x, v}-mode. Then 
the process is said to have a 'short-range' (or short) memory, if there is a finite T* > T 
such that for all t > T* {p 2 ,fi, v}-mode holds. Otherwise, the process is said to have a 
'long-range' (or long) memory. 
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Remark 3.1 According to definition \3. a diffusion mode has a long memory if the 
influence of the old diffusion mode never vanishes, even though the diffusion mode is 
changed, i.e. the particle does not forget its past. In the case of short memory, the par- 
ticle remembers the old mode for some critical time, and then forgets it fully, recognizing 
the new mode. 

Theorem 3.1 Let v > and fi ^ 0. Assume the {Pi, fi,v}- diffusion mode is changed 
at time t = T to the {P2, fi, v}- diffusion mode. Let T* = T / fi and t* =T/(p + v). Then 
the process has a short memory. Moreover, 

(i) {Pi, fi, v}-diffusion mode holds for all < t < t*; 

(ii) {fa, fi, v}-diffusion mode holds for all t > T*; 

(Hi) a mix of both {Pi, fi, v\ and {fa, v}-diffusion modes holds for all t* < t < T* . 

Proof. Let f3(s) = Pi for < s < T and [3(s) = l3i for s > T. Assume v > 0. 
Denote s = fit + vt. So, the {[3\, fi, z/}-diffusion mode holds if fit + vt < T. Let 
< t < t* = T/{fi + v). Then for every r G (0, t) we have fit + vt < {fx + v)t < T. 
This means that the order operator (3{s) in Pfr^ , v y takes the value fix giving (i). If 
t >T/ fi then for all r > 0, fit + vt > T. Hence, f3{s) = fa, obtaining (ii). Now assume 
T/{fi + v) < t < T/fi. Denote tq = {T — fit) /v. Obviously r > 0. It follows from 
{fi + v)t > T dividing by v that t > T/v — tfi/v = r , i.e. < r < t. It is easy to check 
that if < t < tq then fit + vt G {fit, T) C (0, T), giving (3{s) = Pi, while if r < r < t 
then fit + vt G (T, {fi + v)t) C (T, oo), giving P{s) = Pi- Hence, in this case the mix of 
both {Pi, fi, v} and {P2, fi, z/}-diffusion modes is present. ■ 

Theorem 3.2 Let v < and fi+v 7^ 0. Assume the {Pi, fi, v}-diffusion mode is changed 
at time t = T to the {P2, fi,v}- diffusion mode. Lett* =T/fiandT* =T/{fi + v). Then 
the process has a short memory. Moreover, 

{i') {Pi, fi, v}-diffusion mode for all < t < t* ; 

{ii') {P2, fi, v}-diffusion mode holds for all t > T* ; 

{in) a mix of both {Pi, fi, v} and {P 2 , fi, v}-diffusion modes holds for all t* < t < T* . 

Proof. Let v < 0. Assume again P{s) = Pi for < s < T and P{s) = p 2 for s > T. 
As in the previous theorem, denote s = fit + vt. First we notice that if < t < T/fi 
then fit + vt < T, which implies P{s) = Pi, giving {i ). Now let t > T/{fi + v) be any 
number. Then for < r < t we have fit + vt > T, which yields P{s) = P2. So, we get 
{ii ). Now assume T/fi < t < T/{fi + v). Again denote r = (T — fit) /v. Obviously 
tq > 0. It follows from {fi + v)t < T dividing by v < that t > T/v — tji/v = r , i.e. 
< tq < t. It is easy to check that if < r < t then fit + vt G (T, fit) C (T, 00), giving 
P{s) = Pi, while if r < r < t then fit + vt G {{fi + v)t,T) C (0,T), giving P{s) = Pi. 
Hence, in this case the mix of both {Pi, fi, v} and {/^2,/i, ^-diffusion modes is present, 
obtaining {Hi ). ■ 

Corollary 3.1 Let v = and fi 7^ 0. Assume the {Pi, fi,v}- diffusion mode is changed 
at time t = T to the {P2, fi, v}- diffusion mode. Let T* = T/fi. Then the process has a 
short memory. Moreover, 
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(a) for all < t < T* there holds v}- diffusion mode; 

(b) for all t > T* there holds {/^2,/x, v} -diffusion mode. 

Proof. If v = then we have /3(s) = = /3i for t < T/fi and /3(s) = /3 2 for 

t > T//i. ■ 

Corollary 3.2 Let /x = orfi + u = 0. Assume the {Pi, /i, i/}- diffusion mode is changed 
at time t = T to the {P2, fi, u}-diffusion mode. Then the process has the long memory. 

Proof. According to the structure of LH-parallelogram \x = implies v > 0. In this 
case T* = oo. H fi + u = then < and t* = oo. In both cases we a have long memory 
effect. ■ 

Remark 3.2 Notice, that ifu = then there is no intervals of mix of modes. Moreover, 
if v = 0, // = 1, i7ien T* = t* = T. In this sense we say that a process has no memory. 
For all points {(/i, u)} except {/i = 0,0<z/<l} and {v < 0, fi + v = 0}, £/ie operator 
D^j^ v -, has a short memory. The memory is stronger in the region v < and weaker 
in v > 0. On the line \i + v = 1 we have t* = T < T*. The lines /i = 0, v > and 
/i + z/ = identify the long range memory. 



4 The Cauchy problem for variable order differential 
equations 

In this section we study the Cauchy problem for variable order differential equations 
with a piecewise constant order function (3{t) = J2k=o^Pk, where T k is the indicator 
function of [T k , T k+1 ). We assume that the diffusion mode change times T 1; T 2 , ...,T N are 
known, and set T = 0, T N+ i = oo. We assume that the solution of the initial value 
problem for the VOPDE ([TBI) is continuous^] when the diffusion mode changes. 
Thus, the Cauchy problem is formulated in the form 

V^ u} u(t, x) = A(D)u(t, x),t>0,xeM n , (14) 

u(0, x) = <f(x), (15) 

u{T; -0,x) = M(T fe * + 0,x), k = 1,...,N, x G M n , (16) 

where A(D) is a pseudo-differential operator with a continuous symbol A(£), £ G R n , 
and T£ are actual mode change times. It follows from Theorems 13.11 and 13.21 that are 
defined through T k /(fi + v) and T k J 

We note that, since the integration operator order (3 depends on the variable t, a 
variable order analog of the integration operator becomes 

jm f(t) f' (t-r)«-'+->-'/(r) 



3 in the topology of ^ GtP (M n ) (or ^^ GtP )- 
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which we call a variable order integration operator. 

For further purpose we recall the definition of the Mittag-Leffler function [23] in 
the power series form 

oo n 
n=0 v ' 

Obviously, Ep(z) = e z , if j3 = 1. For all (3 > Ep(z) is an entire function of type 1 and 
order 1//3. Note that Ep(—t), t > 0, is completely monotone [23], and has asymptotics 
£/?(-*) = Oft" 1 ), t -> oo. 

Lemma 4.1 Assume < /3* = mino<j<jv{^} and [&/3«] i/ie integer part of k{3*. Let 
v(t) be a function continuous in [0, oo). T/ien /or arbitrary T > and even/ = 1, 2, ... 
i/ie estimate 

max 0<t<T \J^ k v(t)\ < ^ T } ]k . max 0<t<T \v(t)\ (18) 

+ 1J! 

r^*, 0<r<l, 



z/>(r) 



r, r > 1. 



Proof. Let uft) be a function continuous in [0, oo). For k large enough, so that 
/3*k > 2 we have minT(k/3(fit + z/r)) = r(fc/3*). Taking this into account, for all such k 
and for all t e (0, T] we obtain the estimate 

J T(k(3(ixt + VT)) r(A;/3* + l) 

and hence, the estimate in Eq. f[T8"l) . ■ 

Let t cr> j = TjKn + u),j = 1,...,N, be critical points corresponding to the mode 
change times Tj,j = 1,...,N. We accept the conventions t cr>0 = 0,t crN+ i = oo. Let 
Ep(z) be the Mittag-Leffler function with parameter (3 G (0, 1]. Now we introduce the 
symbols which play an important role in the representation of a solution. Let 

Sj(t,0 = E p .{{t-t cr ^A{i)), t > tcrj, j = 0, N, (19) 

and 

fc-i 

M}.(t, £) = Sfcft — t cr ^, £) Y\ Sj(tcr,j+1 — tcr,j, 0) t — tcr,k, k = 1, N. (20) 

3=0 

Further, we define recurrently the symbols 



r(l - A) J (t - r)* r(l - (3i) Jo r^{t - r)ft ' 

P_ift,£) = 0, P ft,0 = 1, ttft.O = ^ $ft + W - r.Ow^^^Odr, 



and if 



^•ft, = / Sift + W - r, t^Dl-^ir, Odr, 

J t c .r ■ 



is defined for t > t cr j and for all j < k — 1, then for t > t cr ^, 



for jfe = 2,...,JV. 

4.0.1 The case v = 0. 

Theorem 4.1 Assume u — and </j G ^>G, P {lR n )- Then the Cauchy problem (Lfy- fTb)) 
has a unique solution u(t,x) G C([0, T], 1 i r Q iP (iR n )) ; T < oo, which is represented in the 
form u(t,x) = S(t, D)ip(x), where S(t,D) is the pseudo- differential operator with the 
symbol 

N 



s(t,o=i' So(t,0 + J2 1 '^ i M ^0 + 

k=l 

S k (t, f ) / S k ^(t^ k + t cr , k ^ - r, t^D^TZ^r, £)dT + 

* tcr,k — l 

j S k {t + t cr , k -r,0 tcrik D}T^n k {r^)dr} , (22) 

wherel' k = T[t crk ,t crk+1 )(t), k — 0, ...,N, are indicator functions of the intervals \t cr ^ k ,t cr ^ k+ i) 
k = 0,...,N; SjHi), j = 0, N, M k (t,£) andK k {t,£), k = 1, N, are defined in / TTPj) . 
UZO) and l[21\) . respectively. 

Proof. It is not hard to verify that 

N N 

^iV^£o}«M = ^2^ k J pk D^ k u(t,x) = u(t,x) - J2^k<t cr , k ,x)+g(t,x), (23) 

fc=0 fc=0 

where 

N , t Pk fc_1 

g(t, X ) = Y, 1 'k V{l+ ^ k) H t^Vt^oyUitcrj+uX). 

Multiplying both sides of equation (fHl) by jfjflx and applying the formula (|23|) . we 
obtain 

N N 

u(t,x)-Y,Z'k jf3k A(D)u(t,x) = Y,?k u (tcr,k,x)-g(t,x). (24) 

fc=0 k=0 

Let t G (0, t C r,i) • Then /3(f£t) = Po and g(t,x) = 0. In this case taking into account the 
initial condition (|15p . we can rewrite equation fl24l) in the form 



■ft 



(t, s) - J 0o A(D)u(t, x) = <p(x), 0<t< t crA . (25) 
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The obtained equation can be solved by using the iteration method. Determine the 
sequence of functions {uo(t, x), u m (t, x)} in the following way. Let Uo(t, x) = <p(x) 
and by iteration 



u, 



,(t, x) = jP°A(D)u m ^i(t, x) + <p(x), m = 1, 2, ... (26) 



We show that this sequence is convergent in the topology of the space C[0,T; \l/(i? n )] 
and its limit is a solution to the Cauchy problem (fl^|) . (fl5|) . Moreover, this solution can 
be represented in the form of functional series 

oo 

u(t, x) = J2 J Pok A k {D)ip(x). (27) 

Indeed, it follows from the iteration process (1261) that 

u m (t, x) = J^ m A m (D)ip(x) + J^ m -^A m -\D)(p{x) + ... + ip(x). (28) 

Now we estimate u m (t,x) applying Lemma [TBI term by term in the right hand side of 
(T281) . Indeed, let iV G IN. Then taking into account the fact that the Fourier transform 
in x commutes with J' 3 ™, we have 

max [0 , T]PN (J^ k A k (DMx)) < \ ^ p N (A<p(x)). 

[k(3 ]\ 

Further, since A(£) is continuous on G there exists a constant > 0, such that 
max^ EsuppKN \A(^)\ < C N , or, by induction max^ suppKN \A k (^)\ < C%. Hence, for every 
N G IV we have 

PN (J^ k A k (DMx)) < Mp "^^ • (29) 
It follows from (ESI) that 

C k N ^ k {T) 



max[ 0> np N (u m (t,x)) < \\<p{x)\\p22 



T((3 k + 1) 



k=0 

<Cy(x)\\ p E Po (C N ^(T)),N = l,2,..., 

where Ep (r) is the Mittag-Lemer function corresponding to /3q. As far as the right hand 
side of the latter does not depend on m, we conclude that u m (t,x) defined in fT28|) is 
convergent. Again making use of Lemma [TBI we have 

p N (u(t,x)-u m (t,x))<\\cp(x)\\ p J2 YWk + lV N = l ^---- (30) 

k=m+l 

The function 

OO L, 

1] 



fe=m+l 
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on the right side of equation f[3"Ul) is the residue in the power series representation of 
the Mittag-Leffier function Ep (r)), and, hence, TZ m (rj) — > 0, when m — > oo for any real 
(or even complex) 77. Consequently, u m (t,x) — > u(t,x) for every N = 1,2, that is in 
the inductive topology of the space C([0, 00), \I/ G)P ). Thus, u(t,x) G C([0, 00), ^g,p) is 
a solution. Moreover, it is readily seen that u(t,x) in (1271) can be represented through 
the pseudo-differential operator S(t,D) with the symbol S (t,^) = Ep (t^°A(^)) in the 
form u(t,x) = Uo(t,x) = So(t, D)ip(x), t G (0,t cr] i). By construction u(t,x) is unique 
and continuous in t. So, lim t ^ tcrl __ u(t,x) = E / 3 (t^ 1 A(D))(p(x) exists in ty G p(lR n ). 
Further we extend u(t,x) to \t CT ^,t cr ^)- Equation (fl4|) in the interval (t cr ,i, t cr ,2) reads 

D^u(t,x) = A(D)u(t,x), t G (t cr ,i,t cr , 2 ). 

Splitting the integration interval (0, t) on the left hand side of the last equation into 
subintervals (0,t crj i) and (t crj i,t), we can rewrite it in the form 

tcrA D^u(t,x) =A(D)u(t,x) + F l (t,x), t G (t cr ,i,t cr , 2 ), (31) 

where 

^■*> = -r(T^)i t^F*' 

Taking into account the fact uo(t, x) = So(t, D)ip{x), it is not hard to see that F\{t, x) = 
IZiit, D)ip{x), where 



Due to the continuity condition (1161) . we have also 

u(t CT) i,x) = u (t crA ,x) = E l3o (t^ 1 A(D))ip(x) = S (t cr ,i,D)(p(x). (32) 

In the general case, assuming that solution are found in the intervals [0, t cr ^), [t C r,fc-i, t cr;k ), 
we have the following inhomogeneous Cauchy problem for the interval [t cr ^ k ,t cr ^ k+ i) 

tcr k D^u{t, x) = A(D)u(t, x) + F k (t, x), t G (t cr , k , tcrj+i), (33) 

u(t crik ,x) = Uk-l(t cr!k ,x), (34) 

where 

It is not hard to verify that F k (t, x) can be represented in the form TZ k (t, D)ip(x) with a 
pseudo-differential operator 1Z k (t,D) whose symbol is given in (pTj) . A unique solution 
to (|33|) . (|34j) can be found applying the fractional Duhamel principle (see [35] ) 

u k (t,x) = S k (t,D)u k ^ 1 (t cr . jk ,x) + 
I Sk (t - (r - t cr , k ),D) tcr k D l -^F k (T, x)dr, 
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k — 1, iV, t crjk < t < t cr ^+i- 
Now taking into account that 

u k-l(tcr,k? = [ J Sj(t cr j+i — tcrjj -D) 

i=o 



C "' fc Skater* - (r - tcr,fc-i), ^) tcnfc _ 1 ^^- 1 F A; _ 1 (T, x)dr 

^cr, fc — 1 

we arrive at fl22l). ■ 



Remark 4.1 Assume in equation [Lfy (3{t) = (3, where (3 is a constant in (0, 1]. Then 
the representation formula (2B) is reduced to 

u(t,x) = E/3(t P A(V))(p(x), 

which coincides with the result of /IT]/ . 

Applying the technique used in the paper [TT] and the duality of the spaces ^G,p(lR n ) 
and ty_Q q {lR n ) we can prove the following theorem. 

Theorem 4.2 Assume v — and <p £ ^_ Gq (lR n ). Then the Cauchy problem (T4\)- (TS\) 
(with '-D' instead of 'D') has a unique weak solution u{t,x) E C([0, T], ^_ G q (M n )), 
T < oo, which is represented in the form 



u(t,x)=l' S (t,-D)i P (x) + J2 I 'k( t ) {M k (t,-D)<p(x)+ 

k=l 

S k (t,-D) [ CrL Sk^it^ + t^^-r^-D) t^Dl-^n^T.-D^^dr + 

J t cr ,k-l 

f S k {t + t cr , k -r,-D) tcrk D 1 -^n k (r,-D)ifi(x)dr}, (35) 

Corollary 4.1 If v = then the fundamental solution of equation [Lfy with the conti- 
nuity conditions in (To)) is represented in the form 

U(t,x)=f (t)-±- [ E Po (t^A(-0)d^+ 

N „ fc-1 

k=l [Zn) jRn j=0 

Ep k ((t-t k f*A(-0) P E^h-rf^Ai-i)) tk _ 1 D 1 T -^n k ^(r,-Odr+ 

f E Pk ((t - r)^A(-O) t k Dl-^TZ k (T, (36) 
Jt k 

where tj = t cr j. Moreover, U(t, x) G ^_ Gq {]R n ) for every fixed t > 0. 
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4.1 The case -1 < v < 1 



The solution u(t,x) = S(t, D)tp(x) obtained in Theorem 14.11 in the case v = has 
the structure u(t,x) = D)(p(x) + ^ 2 (t, D)ip(x), where ^i(t, D) and ^ 2 (t,D) are 
operators with symbols 

JV 

* ! (t , o = i' s (t , o + z'k (t) M k (t , o , 

fc=i 

and 

* a (U) = 

^ P^cr,k 
fc = l "'icr.fc-l 

/ SfcCt + tcifc-T.O ^^^^(r^)^} (37) 

*^ ^cr,k 

The term v(t,x) = ty(t, D)tp(x) reflects an effect of diffusion modes, while the term 
w(t,x) = ^ 2 (t, D)cp(x) is connected with a memory of the past. We note that this 
structure remains correct in the general case v G (0, 1] also, however, the symbols of 
solution operators are further restructured depending on the intervals of mixture of (two 
or more) modes. The theorems formulated below concern some intervals free of mixed 
modes. 

Theorem 4.3 Assume /i^O, + and tp G ^^(iR™). Then there exists a number 

T* > and pseudo-differential operators V*(D) and lZ*(t,D) with continuous symbols, 
such that for t > T* the solution of the Cauchy problem [Lfy- fTb]) coincides with the 
solution of the Cauchy problem 

T ,D^ N u(t, x) = A{D)u{t, x) + /*(*, x), t > T*, x G R n , (38) 

u(T*,x) = (p*(x), x G JR n . (39) 
where f*(t,x) = K*{t, D)tp(x) andip*(x) = V*{D)ip(x). 

Proof. Assume v > 0. Then as it follows from Theorem 13.11 that the actual mode 
changes occur at times T* = Tj/ \i and t* — Tj/(fi + is), j = 1, ...,N, if diffusion modes 
change at times Tj,j = 1,..,N. Obviously, t\ < ... < t* N and T x * < ... < T£ if T x < 
... < Ijv- The order function /3(fit + vt) under the integral in D^Pj, takes the value 
/3jv for all t > and r > 0. Hence, the variable order operator on the left side of 
ffl4l becomes D% N if t > Tjy. Analogously it follows from Theorem 13.21 that if v < 
then /3(fit + vt) takes the value (3n for all t > t* N and r > 0. Thus, if v ^ then 
for all t > T* = max{T^,t* N } and < r < t we have (3 (fit + vt) = (3n- Similar to 
the case v = 0, splitting the interval (0, t), t > T*, into subintervals, we can represent 
the equation f fT4"l) in the form f l38i) . Further, from the continuity condition ( fl6l) we 
have u(T*,x) = lim 4 ^T*-o v(t, x), where v(t,x) is a solution to the Cauchy problem 
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for fractional order pseudo-differential equations in sub-intervals of the interval [0, T^) 
constructed by continuation. Therefore there exists an operator S*(t,D), such that 
v(t,x) = S*(t,D)<p(x). Denote V*(D) = S*(T*,D). Then u(T*,x) = V*(D)ip(x). This 
means that for t > T* solutions of problems (fT4|) - f|T6|) and (!38l) . (l39l) coincide. If is = 0, 
then the statement follows from Theorem 14.11 ■ 

Theorem 4.4 Assume (p G \E f G i p(iR n ). Then there exists a number t* > 0, such that for 
< t < t* the solution of the Cauchy problem [Lfy- fTb}) coincides with the solution of 
the Cauchy problem 

D^u{t, x) = A{D)u{t, x), t > 0, x G FT, (40) 

u(0,x) = (p(x), x G FL n . (41) 

Proof. It follows from Theorems 13. II and 13.21 that the order function (3(fit + isr) under 
the integral in takes the value ft for all t < t* = mm{t*,T*} and < r < t. 

Hence, the variable order operator in (fl4l) becomes Dp if < t < t*. The order ft 
(or diffusion mode {ft, fi, is}) has no influence in this interval. For t > t* two diffusion 
modes {ft), /•*, v} and {ft, /i, is} are present. If t < min{t\, T 2 * } then for all r > we have 
fit + 1st < T 2 . That is, there is no influence of the mode {ft, fi, is} if t < min^^T^}. 
In the same manner the other values of f3 have no influence in the interval < t < t* . 
This means that for < t < t* solutions of problems f|T4l) - f|T6l) and (J1D]1 . (j4T]l coincide. ■ 



5 Some properties of the fundamental solution 

In this section we study some basic properties of a solution u(t,x) to the problem (fUj) - 
(jl~6l) . Namely, we show that u(t, x) is a density function under a rather general conditions 
on the pseudo-differential operator A(D). We also study the MSD(t) of the correspond- 
ing process near the initial time, which is important, in particular, in cell biology. 

Theorem 5.1 Assume A(£) is a continuous symbolwith negative values for £ ^ 0, A{— £) = 
A(£) and A{fS) = 0. Then the Fourier transform U(t,£) = F[U)(t, £) of the fundamental 
solution U(t,x) to the problem [Lfy- fTb}) satisfies the following conditions: 

1. U(t,£) is continuous in £ for every fixed t > 0; 

2. U(t, 0) = 1 for all t > 0; 

3. U(t,£) is positive definite for every fixed t > 0. 

Proof. First, let is = 0. Then Corollary 14.11 and the symmetry A(— £) = A(£) imply 
that 

U(t,0=l' (t)E^A(0)+ 

£ J fc w {E, k ((t - t k y*A(o) n^((*i+i - *i)M0)+ 

fc=i i=o 
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f Ep k ((t-rf*A(0) th D^n k (r,Odr}, (42) 

where Ep(z) is the Mittag-LefHer function and Tlk{t,£) is defined in (j2"Tl) . Taking into 
account continuity of Ep(z), we conclude that U(t,£) is continuous for every fixed t > 0. 
Further, it follows from the definition of lZk(t,£) that 7£(£, 0) = 0. This implies U (t, 0) = 
SfcLo-^fc(*)> smce ^a(O) = 1- Finally, as is known [231130], the function E@{— Xt 13 ), < 
/3 < 1, is completely monotone for all t > if A > 0, and therefore Ep{—\t^) > and 
dE p {-\tP)/dt < for all t > 0. The latter implies ft fc (t,£) > 0, k = 1, JV. It follows 
from this fact together with < and positiveness of Ep(t@ A(£)) for all t > that 
U(t,£) > for every fixed t > and (6XCC where is an arbitrary compact. 
Now it is easy to verify positive definiteness of U(t,£) for each fixed t > 0. The idea of 
the proof in the general case v G (—1,1] is preserved, since the general structure of the 
representation formula for the fundamental solution remains unchanged. ■ 

Corollary 5.1 Under the assumption of Theorem \5.1\ the fundamental solution U(t,x) 
to the problem [Lfy-fHfy is a probability density function for each fixed t e (0, oo). 

The proof of this statement immediately follows from the Bochner-Khinchin theorem 
(see, e.g., [3]). 

Thus there exists a stochastic process X t with a density function pt(x) = U(t,x) for 
every fixed t > with po(x) = S(x). Denote by \i t = E[X t ] the expectation of a random 
variable X t (t is fixed) and MSD(t) = E[\X t - fi t \ 2 ]. 

Now assume that the pseudo-differential operator on the right hand side of ( 1T4"1) is a 
negative definite second order homogeneous elliptic operator, that is the symbol of the 
operator A(D) has the form A(£) = a ij£,i£,j- The matrix A = (a^) is symmetric and 
negative definite: V T AV < — C\V\ 2 , C > 0, where V is an n-dimensional vector, V T is 
its transpose. By Tr(A) we denote the trace of A: Tr(A) = X] a ii- ^ is not hard to 
verify that in this case U(t, —x) = U(t,x) and \i t = 0. Hence, 

MSD(U;t)= [ \x\ 2 U(t,x)dx. (43) 

Theorem 5.2 Assume A(D) in [Lfy is a second order homogeneous negative definite 
elliptic operator. Then there exists t* > such that for t < t* the function MSD(U ; t), 
where U(t, x) is the fundamental solution of the Cauchy problem [Lfy- fTb)) , is represented 
in the form 

Tr( A) 

MSD(U; t) = jv^-f)^ < t < t\ (44) 

Proof. The proof is an implication of Theorem 14.41 and the fact that MSD(U; t) for a 
solution of the Cauchy problem (140p . (14ip can be represented in the form MSD(U;t) = 
(-Ag)#(U)|g=oi 
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Theorem 5.3 Under the condition of Theorem \5."A for MSD(U ; t), where U is the fun- 
damental solution of the Cauchy problem [Lfy- fTb}) , the asymptotic relation MSD(U; t) = 
0(t PN ), t oo, holds. 

Proof. It follows from for t > that U(i,£) = Qi(t,£) + Q 2 {t,£), where 

JV-l 
j=0 

/ E^dtN - r)^-M(O) t^Dl-^n^T^dT }, 

J t AT— 1 



and 



Q 2 (t,0= [ Ep N ((t-rf»A(0) tN Dl-P»K N (T,OdT. 
Jt N 



It is not hard to verify that Qi(t,£) dominates in terms of asymptotics for large t. We 
observe the same after applying — Ag as well. Elementary, but tedious calculations show 
that (-At)Qx(t,g) = 0(t^ N ), t -> oo. ■ 

Corollary 5.2 Let (3{t) = (3, where (3 is a constant in (0, 1]. Then 

MSD{U]t) = WVi) t ^ t>0 - 

Remark 5.1 A natural generalization of the model is to allow changing of random dif- 
fusion modes (3k at random times with appropriate distributions respectively. The 
questions on an asymptotic behaviour of U(t, x) for large times, which shows how heavy 
is the tail of the distribution, as well as an asymptotic behaviour of MSD(U;t), t — > oo, 
which tells about the nature of the corresponding process, are important. We will discuss 
these challenging questions in a separate paper. 
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